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We use Chiral Perturbation Theory to compute the nucleon mass-shift due to finite volume and 
5h I temperature efi'ects. Our results are valid up to next-to-leading order in the "e-regime" (mL ~ 

^ , m/3 <^ 1) as well as in the "p-regime" {mL ~ m/3 S> 1). Based on the two leading orders, we 

discuss the convergence of the expansion as a function of the lattice size and quark masses. This 
result can be used to extrapolate lattice results obtained from lattice sizes smaller than the pion 
cloud, avoiding the numerical simulation of physics under theoretical control. An extraction of the 
low-energy coefficient C3 of the chiral Lagrangean from lattice simulations at small volumes and a 
"magic" ratio P — 1.22262L might be possible. 

m ■ 
> 

CTN . Lattice QCD simulations are necessarily performed in finite boxes. Finite-size eff^ects are controlled by the parameter 
mL, where L is the lattice size and m the mass of the lightest particle, in QCD, the pion. Physical results can be 
obtained in the limit mL 1 . As the pion masses achieved in simulations approach the physical value it becomes 
' harder to fulfill this condition. However, most of the configurations in large lattices describe pions traveling at large 
' distances of the order of L. Since the physics of these soft-pions is strongly constrained by chiral symmetry, strong 
[ theoretical control over them makes their numerical simulation unnecessary. One can thus obtain physical results by 
simulating in smaller lattices and using Chiral Perturbation Theory (xPT) or some other relevant effective theory to 
include the soft-pion physics cut off by the box size and extrapolate the results to the infinite volume limit. Another 
way to describe the same procedure gives added insight: The low-energy physics in the infinite and finite volume are 
(jj \ described by the same effective theory with the same low-energy constants, since the values of these constants 
' encapsulate short-distance physics that is not modified by finite-volume effects. The comparison of finite volume 
] lattice results with the effective theory prediction allows one therefore to determine the value of some of the low 
energy constants. Those, in turn, can be used to determine physical observables in the infinite-volume limit. 
This general procedure has been carried out in the regime mL 3> 1, where standard xPT can be applied, to a 
variety of one nucleon observables, see e.g. 0. However, it is for mL ^ 1 (in the so-called e-regime that the 
programme described above is fully realized. For such small boxes, most of the pion cloud surrounding a baryon is 
excluded, and we are left with a bare nucleon. There are some modifications to the usual xPT power counting in 
this regime. The first and obvious one is that the momenta are quantized in units of 2tt/ L. More importantly, the 
pion zero mode fluctuations are not suppressed, become non-perturbative and need to be treated exactly 2]. They 
reduce the value of the chiral condensate and make the chiral condensate disappear altogether in the chiral limit. 
This is to be expected since there is no chiral symmetry breaking at finite volumes. Recently, the e-regime in the 
meson sector and its relevance to lattice QCD have been assessed in a number of papers W\. In the present work, we 
extend the idea to the one-baryon sector. Convergence in the baryonic sector is typically worse than in the mesonic 
sector, as it receives contributions at every order in p/{Anf), unlike the meson sector case where the expansion 
parameter is (p/(47r/))^. We address this issue by comparing the sizes of leading and next-to-leading order 
contributions in a calculation of the nucleon mass. 
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We consider one nucleon in a small box of size f3 x for 27r/(47r/) ^ f3,L ^ 27r/m, the "e-regime". L is the size of 
the spatial directions, P the temporal extend of the box, namely the inverse temperature. In this regime, xPT is 
valid, except that the relative counting between p and m is changed. Instead of the usual counting 
1/L, 1/ (3 ^ m ^ p (p-counting), we use 1//3 ~ ruq ~ ^/m ^ e, hence the name "e-regime" For small boxes, 
the first non-zero pion mode has a momentum p = 2tt/ L > A, so we include the A(1232) as explicit degree of 
freedom, counting, in the e-regime, A ^ m ~ e^. 

THE e EXPANSION IN THE BARYON SECTOR 

Low-energy properties {Q ^ 1//3) of the system are described by the effective Euclidean Lagrangean 



+4(C2 - ^^)Al + Ac^A^A^ - (C4 + -^)2ie^=^AA,a'^ + ■ ■ -W, 

Ca - -At'^po + A - ^)A'^ + 5iVAAt'^K^iV + H.c.) -f • • • , (1) 

where we list only the terms pertinent to our calculation. The pion decay constant is / = 92.4 MeV, £,l,£,r are 

SU{2) matrices parameterizing the chiral SUl{2) x SUr{2) group and M — diag(TOq,mg) is the quark mass matrix 

in the isospin limit (the precise conventions used can be found in Appendix A). The values of the other low-energy 

constants will be given when we discuss our results. The Goldstone bosons belong to the coset space 

[SUl{2) X SUr{2)]/SUl+r{2), and we are free to choose an arbitrary member to be the representative of each coset 

("fix the gauge"). Instead of the usual choice £,l = £,r = ^ = e^T", we use the choice made in background field 

calculations 

Cr = ule-'V, (2) 
where mq is a space-time independent field and -k^x) does not contain zero-modes: 

^{x)= e'^*+'^-: (3) 

The rationale to separate zero- and non-zero modes is that, as we will see below, the zero modes obey a different 

power counting from the non-zero ones at small volumes where chiral symmetry is partially restored. 

The background field uq only appears in those terms of the action which include quark masses. This can be easily 

seen by noticing that a non-trivial background Uq corresponds to a chiral rotation of the vacuum one expands 

around. In the absence of quark masses, all such vacua are equivalent, so the physics of the Goldstone bosons is the 

same. The terms which do however depend on the quark masses are in the isospin limit with 

7^]Tr(^) = iTr(A-F At): 

- mqBfn\ Tr(u^e^ ) - AmqBciN'^neTi{ule'^ )N. (4) 
At leading order, rn^ = 2Bmq is the pion mass in the infinite volume limit. 

We can now estimate the different terms of the Lagrangean. The typical fluctuations of the non-zero modes 7r(a;) are 
of the order 7r(x) ~ e since, for larger values of 7r(a;), the kinetic term is much larger than one and suppresses their 
contribution to the path integral (we can estimate the size of the kinetic term as 1/e^ coming from the volume 
integral, e^ from the two derivatives and e^ from the pion fields). A similar argument implies TV ^ e^^^. However, as 
observed by Gasser and Leutwyler 0, the zero- mode uo is of order e*'. We can conclude that by noticing that the 
coefficient of the first term of Eq.l^J is of order ~ m^f^(3L^ ~ e°. Because the zero-mode is not suppressed, it has to 
be treated exactly. This is related to the restoration of chiral symmetry at finite temperatures and volumes. In 



3 



small boxes the zero-mode fluctuates over the whole group manifold, in contradistinction to the infinite volume limit 
in which the zero-mode makes only small fluctuations around a preferred vacuum direction. As shown in |2j the 
integration over the zero-mode can be performed as follows. The part of the partition function which contains Uq 
can be written as 

Z[N,A] = J [Duo] exp J d^x {^^^ + 2m'ciN\x)N{x)^ (^eTT{ul) ~ + j^UeTiiu^iTr^T^ ) + 

« J [Duo] e^^^T'("o) 1 + TZeTiiul) (^2m\, J d^x N (^1 - - J d'^x 



X{s) cxp[- 



2X{s) 



where we dropped higher orders in the pion-fields, s — Bf^mqPL 
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P(3L^/2 and 
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(5) 



(6) 



with Ii (x) a modified Bessel function. The integration over the zero- mode performed above renormalizes the 
nucleon mass (adding a term proportional to ci of order to it) and the pion mass (by a term of order e"), as well 
as the non-derivative couplings: 



M 



M - Am^ci 



m^g = 2mqB 



2 X'((3L^m^P/2) 



2X(/3L3m2/2/2)' 
X'{(3L^m^Pl2) 
2X{(3L^w?p/2y 



(7) 



The effective pion mass is shown in Fig.^ In the limit s oo, one retrieves with X'{oo)/{2X{(X))) = 1, the 
well-known infinite-volume results. Notice that the shift of the nucleon and pion masses due to the zero modes 




FIG. 1: The effective pion mass rricff as function of PL^m? p /2, Eq. (|7J. 
comes just from quenching the chiral condensate in the finite volume: 

(0 |gg|o)/3.L _ x'{s) 

{{)[qq[Q) 2X{sy 

The total partition function of the system is finally to the order considered 

Z = j [DN] [DA] e^^' Z[iV, A], 
with S' the part of the action which is independent of zero-modes. 
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FIG. 2: Leading-order contributions to the nucleon mass. 



NUCLEON MASS 



The shift in the nucleon mass due to finite volume effects is given at leading order [©(e'^)] by the two one-loop 
diagrams of Fig. |21 We find for the first diagram: 



ifi|A(A-0,meff). (10) 



.0 .V 



The second diagram is 



A 2 ■ f27Tn\2 



3/' „'^,^+2^+^A(2^Iia)2 + (2^)2+, 



eff 



^' -^^A(A,meff). (11) 

Because mL, m/3, AL and A/3 are all of order e in the e expansion, the m and A contribution to these graphs are of 
order O(e^), so that m (and with that of course TOcff) and A can be dropped from the expressions above. 
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FIG. 3: Next-to-leading order diagrams for the nucleon mass. The square vertex represents a vertex suppressed by e, the cross 
a kinetic energy insertion, and the circle the zero-mode mass contribution. The dashed, full and double lines represent a pion, 
nucleon and A propagator, respectively. 

The first truly specific feature of the "e-regime" appears at order e^, because the nucleon mass receives contributions 
from the zero modes computed above in Eq. ((TJ in addition to the graphs shown in Fig. |3| The first graph leads to 

^(4)/^ ^^ ^ 3C2 + 3C3 1 ^ J_ V (¥)^ 

"^ ^^ /2 ;3L3 2^^(2^)2 + (2^)2+„,2^ /2 ^3^3 (2.^)2 + (2^)2 + ^2^ 

^^"^ + 3^-^-'™.ff)-§D(meff), (12) 



T^t,Q o^r 

8M 



with C2 — C2 — 44y. The second and fifth graph vanishes. The third one gives 



M(4)(/3 l) = ^ — V 



o 2 1 / ■ \ /27m\4 
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1{A = 0,meff). 



The fourth graph is the non-perturbative contribution computed before in Eq. ((TJ as 



2 X'(mV^/3£V2) 



(13) 



(14) 



and the last one contributes as 



Mf^(/3,L) 



2.9^ 



2 



1 



3M/2 



E 



(2^)4 



^^B(A,meff). 



27rno 



+ iA 



M 

3M/2 



(15) 



The functions A,B,C and D are calculated in Appendix B. We reduced them to rapidly converging sums for 
non-zero values of rrics and A, but no analytic form is available. A Mathematica notebook computing these 
functions is available from the authors' website ^. In the e expansion, the contributions coming from the finite values 
of m and A appear only at order in the loop diagrams, so we should take for these to = A = at the order e** we 
are working. In this case, a simple form for the nucleon mass-shift is available: 



,,.(3+4) _ 1 fS9\ . ^9'na\ / Tif3/L) \ 3c-2 
/2i3 ^ 8 3 ; V ML J pL^ 



-2to2ci 



X'{m'^p[3L^/2) 



- 2 



where 



1 1 



1.675 



2vrj 



(16) 
(17) 



with j — \j\. The function t{x) is plotted in Fig.^ 
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FIG. 4: t[(3/L) as a function of the box asymmetry P/L. 

For not-so-small boxes satisfying mL ^ /3A ^ 1, the p expansion applies. The two leading orders in the expansion of 
the nucleon mass in the p-regime are very similar to the ones in the e-regime. The differences are: i) leading and 
next-to-leading order are switched as the quark mass insertion proportional to ci is the leading {p^) order 
contribution while the diagrams in Fig. [21 are the next-to-leading {p^) contribution (terms proportional to C2,C3 are 
even higher, namely p**); and ii) the non-zero value of m and A should be kept in the diagrams. For this reason, if 
we keep the pion mass in our calculations, which in the e-regime is a sub-leading (e^) effect, our expressions will be 
valid in both regimes and, in particular, in the intermediate region (L, /?) = 1/m. This way, they also include some, 
but not all 0{e^) pieces ^ . Furthermore, since the pion mass in the e-regime has a correction of order e*^ coming from 
the integration over the zero-mode (Eq.jTJ), we use TOeff in the one- loop diagrams. 



^ http://nta0.1bl.gov/~'bedaque/index.html or http://ph.tum.de/~hgric 

^ Our results are also valid in the limit /3 — > oo as long as this limit is taken at fixed m. If is kept fixed instead with s = m^/3/^L^/2 ~ 
1, the mass term does not prevent the (n = 0, no ^ 0) modes to have large fluctuations and they become non-perturbative. This is the 
5 regime discussed first in Ref. 0. 
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EUCLIDEAN TIME AND THE CORRECT ANALYTIC CONTINUATION 



There is a subtlety in computing the nucleon mass using the combination of the finite-temperature imaginary time 
and heavy baryon formahsms we used. To see that, consider the derivation of the heavy baryon Lagrangean. One 
starts from the relativistic nucleon field ip and performs a field redefinition which reads in Euclidean space 

^(r, f) = e~*^^(7V(r, r) + H{t, f)), (18) 

where M is the heavy nucleon mass and TV and H are the nucleon and (anti)-nucleon fields satisfying 

7°A^ = N^^^H = —H. An "on-shell" ip field has a fast variation with time (9oV' ^ M), while an "on-shell" N 

satisfies doN ~ 0. The Lagrangean in terms of these new variables is 

^5(^070 + Af)V' N\da + ■ ■ ■)N + H^Oq - 2M + ■ ■ ■)H + ■ ■ ■ . (19) 

The "heavy" field H can then be integrated out and we are left with the usual heavy baryon Lagrangean. Notice 
that the anti-periodic boundary condition in the time direction for the relativistic field implies a different boundary 
condition for the heavy-nuclcon field 

V;(/3, r) = -i^{0, r) ^ N{f3, r) = -e'^'^^N{Q, r). (20) 

Therefore, the field A'^ has the Fourier decomposition 

7V(r,f) =^e-*(^^^^+'^'^)"iV(no,f) . (21) 



no 



The correlators of the field N are defined only at shifted values of (imaginary) frequency, namely at 
oj = Tr{2n + l)/P + iM. 

Consider now, as an example, the computation of the first diagram in Fig. |21 For simplicity, we consider the infinite 
spatial volume limit. As shown in Appendix C, the sum over n can be performed resulting up to constants in 



^+UJ + lA (2|«)2+^2 



13 

dk^ 1 P . f dk^ fc2 I 



nY ko+uj + iAk:^+Lol J (27r)3 (w + iAf + ujI e'3('^-»'^) 
dk^ P uj + iA. 1 



(27r)3 (w -I- i A)2 -I- ojI - 1 " 

where uj is the external energy and Lof. = k"^ . We now substitute lu from above into the second term. 



(22) 



leading to the correct statistics for fermionic ensembles. Physically, that we neglect these fluctuations just mirrors 
the fact that finite-temperature fluctuations of heavy particles are much smaller than those of light ones for 
temperatures (SM 3> 1 at which the heavy-baryon formalism applies. Therefore, we drop this term and arrive at 

/■ dfc4 1 P r dk^ P uj + iA 1 

^^'^^ " J (27r)4 ko + LU + iAk^o+ujl^ J (2^)3 {lj + iA^ + luI ojk e^-- - 1 ' ^ ' 

The nucleon propagator at any value of the external energy (including real values) can be obtained from the 
expression above by analytically continuing in w. In particular, the value determining the mass is obtained for lj ~ Q. 
Clearly, this procedure seems arbitrary for two reasons. First, it seems to depend on the order between setting oj to 
w = 7r(2n + !)/(} + iM and analytically continuing to w = 0. Second, the knowledge of the propagator at discrete 
values of the frequency is not, in general, enough to determine the propagator on the whole complex plane. One 
could, for instance, have maintained e"*'''^ instead of substituting it by —e^^ and usinge"^^*" — l. Fortunately, 
the analytic continuation is unique for functions vanishing at infinity at least like l/|w| pj, as in the case at hand. 
Still, to confirm that we have picked the correct analytic continuation, we repeat this calculation in Appendix C 
without using the heavy baryon formalism in another method to compute finite-temperature corrections which does 
not require an analytic continuation to the real axis, namely the "real time formalism" . 
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NUMERICAL EXAMPLES AND DISCUSSION 

We now present some numerical examples in order to explore the convergence of the e-expansion and to discuss how 
it can be used in the one baryon sector. The leading order result depends on two low energy constants qa = 1.267 
and ^ATA, as well as on the masses and mass splittings m, M and A, whose experimental values are reasonably well 
known. At next-to-leading order, the constants ci,C2 and C3 appear. They are determined experimentally through 
the analysis of pion-nucleon scattering. As they are most sensitive to the isoscalar part of the amplitude, where 
different phase shift analyses disagree, large uncertainties exist in their determination. 

When considering very low energy observables, as we are here, the inclusion of the A as explicit degree of freedom is 
optional. Let us first discuss the case where the A is included. In this case, a determination of the low energy 
constants was made by comparing calculations to the pion-nucleon phase shift data 0|. In this work, different fits 
were discussed using two different phase shift analyses, and also including information about the cr-term. The values 
of ci and C3 are more stable among different fits, while ci varies much more. 

Eq.lUni) shows however that for a certain value of the ratio /3/L « 1.22262, t(1. 22262) = 1/1.22262 and the 
dependence on ci disappears. Since ci and gwA are much better determined, one might use the mass-shifts measured 
around this ratio on the lattice to determine C3. The C2-contribution is generically negligible for /3/L ~ [1 . . . 1.7]. 
In Fig. [SI we present results using the parameter set (fit ^ of Table 4 in Q) : 

5ArA = 1.00 ±0.08, 

ci = -0.35±0.09/GeV, 

C2 = -1.49±0.67/GeV, 

C3 = 0.93±0.87/GeV. (25) 

The errors quoted come from the fit and under-estimate the uncertainty in the constants from higher-order 
corrections. In the left panel of Fig. we show the leading contribution and its next-to-leading order correction to 
the mass shift computed both using Ea. (|16|l and the full formula with finite We// and A. The expansion in 
m/(47r/), A/(47r/) seems to converge for the low-energy constants ci in the given range, except for those close to the 
upper limit and for L smaller than about 1.5 fm, where one approaches the breakdown scale: L/{2-k) « l/(47r/). 
The right panel displays the total mass shift up to second order for different values of C3 between the minimum and 
maximum values suggested by Ea. (|25|l . We notice that the cancellation of the C2-contribution for (3 = 1.22262L 
works very well even for non-zero pion masses, with its contribution to the mass shift never exceeding 5 MeV even 
for m = 300 MeV. 




1 1.5 2 2.5 3 3.5 4 1 1.5 2 2.5 3 3.5 4 

L [fin] L [fin] 



FIG. 5: Left: Finite-volume mass-shift of the nucleon in the theory with explicit A degrees of freedom in MeV as function of 
L [fm] with the central values of the parameter set in Eg. 12511 . Leading order with full expansion (blue solid line); and using 
m = A = Eg. 1161 (green dash-dotted). Next-to- leading order correction in the full expansion (red dashed); and from Eg. ll6ll 
(green dotted). The gray zone shows the variation of the mass shift as C3 varies in the range given in Eg. 12511 . with the upper 
limit corresponding to C3 = 1.8/GeV. Right: Total mass-shift at leading order (blue solid line) and at leading -|- next-to-leading 
order for C3 = 0.06/GeV, 0.93/GeV and 1.8/GeV from top to bottom (red dashed). The parameters (3 = 1.22262L, m = 100 
MeV and A = 294 MeV are the same for both figures. 

In a effective theory without explicit A, its large role in the pion-nucleon interaction is absorbed by the couplings C2 
and C3. In fact, a simple tree level model of the A contribution gives a contribution of 

C2 = — C3 = .g^A/(2(A^ — rn?)) ~ 4/GeV. These values are somewhat larger than what is expected from naive 
dimensional analysis arguments and puts the quark mass expansion in check. The values suggested by different fits 
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roughly agree with the A saturation estimate. We show in Fig. (BJ as an example, the mass shift computed 
with the central values of the parameter set 



Cl = -0.81±0.15/GeV, 
C2 = 2.99±0.77/GeV, 
C3 = -4.70±0.95/GeV 



(26) 



advocated in 0, as well as for a somewhat smaller value of C3 — — 3.4/GeV found e.g. in the partial wave analysis of 
nucleon-nucleon scattering Q. The convergence is obviously poor in either case. While the contributions from ci 
and C3 can be made small for certain ratios this cancellation depends sensitively on the particular values 
chosen for ci and C3 and is hence less useful for lattice determinations. At the ratio (3/L = 1.22262, the 
C2-contribution disappears as before, but it is already negligible at /3 = L. Indeed, a plot of the next-to-leading order 
correction with (3/L — 1.22262 differs from Fig.Elat most at the 5%-level. 
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FIG. 6: Example of the finite-volume mass-shift of the nucleon at leading order (solid line) and the next-to-leading order 
correction (dashed lines) in MeV as function of L [fm] in the theory without explicit A degrees of freedom. The lower dashed 
curve uses the central values of the parameter set in Eg. 12611 . the lower one changes only C3 = — 3.4/GeV following 0|. The 
pion mass is m = 100 MeV, and /3 = L. Dot-dashed (dotted): LO (NLO) for m = A = 0. 



One might wonder why the expansion in the e-regime is so sensitive on the low-energy constants and converges 
badly in the case without A, while it is well behaved in the infinite volume limit. This seems to arise because ci and 
C3 appear only at third order in the p-expansion and are, consequently, poorly determined. In the e-regime, they 
contribute already at second order, so that the uncertainty in their values is enhanced. 



CONCLUSIONS 



We have computed the nucleon mass in a finite box of size j3 x satisfying 47r/ ^ 27r//3, 271 /L^ m (e-regime). 
Taking the value of the low energy constants suggested by experiment, we find that the expansion seems to converge 
for the values of the low energy constants allowed by fits made to pion-nucleon scattering data, if the A(1232) is 
taken into account as explicit degree of freedom. We notice that a particular shape of the box (/3 « 1.22262L) 
eliminates the dependence on one of the low energy constants (02) and suggests determining the value of the most 
poorly known one (C3) by a fit to the nucleon mass with a few different box sizes. We also discussed a subtle point 
involving the combined use of the imaginary time and heavy baryon formalisms. 

We close with the remark that it may seem strange to compute pion loops in small boxes of sizes « 1 fm^, since the 
momentum of the first non-zero mode p « 1.2 GeV is well above the range of validity of Chiral Perturbation Theory. 
An alternative to our procedure that is not subject to this criticism would be to integrate out all modes but the 
zero-mode and obtain an effective theory which is valid only for zero-momentum observables like the mass computed 
here. We point out, however, that in order to connect the low-energy constants of this new theory with the 
parameters of the original chiral Lagrangean, one has to perform a matching calculation that is equivalent to the 
calculation presented in this paper. 
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APPENDIX A: CONVENTIONS 



We collect here the definitions used in the construction of the chiral Lagrangean. Elements of SULip^) x SUr{2) are 
parameterized by {^l,^r)- N is a. spin and isospin doublet with its indices not shown explicitly. A*^ is a 
spin-isospin 3/2-field, where the vector and isovector indices jA are explicitly shown, besides the implicit spin and 
isospin indices. In addition, it satisfies cr^A*'^ = r^A*^ = so only 4 spin and 4 isospin entries are independent, as 
expected for (iso)spin 3/2 objects. The chiral transformation rules are 

N hN, 

A'^ ^ ^Tr(/i-V^/iT^)/iA'^, (27) 



QAB 

where is an orthogonal matrix which is determined by h via h~^T^h = O^^t^ , where h in turn depends on 
Tr{x), L and R. In addition, we define some objects with simple chiral transformation rules 

= l{^lid^^R + ild^^L)^hV^h-'+hd^h-\ 

A = l{^Rd^^R-dd^^L)^hA^h-\ 

D^N = d^N + V^N ^ hD^N, 
D^A'^ = d^A'^ + VpA'^ + ie'^^^Tr(T^V^)A*^ ^ O'^^hD^A'^, 

= Tr(r^yl^)-O^X' 
w^^ = Trir^DM - (28) 



APPENDIX B: CALCULATION OF SUMS 

In this appendix, we drop for clarity the subscript of the efi^ective pion mass nies and denote it by m. The LO 
diagrams contain then the sum 



A(A,m) = — 



/3L3 2^ uj+^+ lA )2 + (2^)2 + „2 



= Ao{A,m)+Ap{A,m), (29) 

where co = 2Tr{k + l/2)//3 is a discrete external energy and k an integer. The ultraviolet divergence of A is identical 
to the one in the infinite-volume diagram and cancels in the difference between finite and infinite volume masses. 

For this cancellation to occur, it is important to use the same regulator in both calculations. In practice, we should 
therefore define the sum above using dimensional regularization. Aq and A^ are the temperature independent and 
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the finite temperature parts: 



lAk^ + (2^)2 



A^(A,m) - ^2 _ (A _ jc^)2 ^„ e/3^„ - 1 " ^ ^2 - (A - iuj)^ e'3(A-»-) _ i ' ^^^^ 

where — ^2™. ^ rn? ^ We used the formula 



'P / ^ y 2tT ~ '^'e*'^^ - j^jllowcrplanc + «R.es(-— gj— |uppcrplano, 



(31) 



which holds if /(z) has no poles on the real axis. We substitute l/{e^^^ ''^^ — 1) by —l/(^e^(^+^^) + 1) sa as in 
Eq. H23|l. The zero-temperature part can be computed with the help of the relation 

LsZ. (2^)2 + ^2 - J {2tt)^ + ^ ^ 2tt 



(fik fc2™ f d^k k 



(27r)3 A:2 _|_ ^2 ' y (2,^)3 ^ fc2 + ^2 



iLk-j 



(Pk fc2™ ^ 1 /"^^ fc(2m+l) sin(jfci) 

^ 9Tr2 r, /„ ; 



where m is a positive integer. Applying this relation to Aq yields 



(27r)3 k^+x^ ^ Jq k"^ + x^ j 

(fk fc2"^ (-x^)" ^ 



Mo = Ao - A(;3 ^ oo, L ^ oo) = ^ ^ i 1 dko'^^e-^^V^ . (33) 



Asymptotically, the sum over j converges because 



,5/2 



SOL.O.A) ^3^ + -- (34) 

These asymptotic forms are also useful in the numerical evaluation of the sum over j. 

We also need Aq evaluated at A = 0. We can obtain this limit noticing that the integral defining g{jL, m, A) is 
nearly infrared divergent when A — s- 0, and hence is dominated by small values of fco- The A — > limit of 
g{jL,m,A) is given by 

A-o 1 f°° „ fc|_+m2 7rm2 
'4 



g(jL,m,A) 



Using this result, 



^o(A = 0,m) = ^;^^. (36) 
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The expression above agrees with that of Ref. 0] . The Hmit m ^ is found by noticing that for small values of ; 
the sum is dominated by the large j terms, which in turn can be approximated by an integral 



(5Ao(A = 0,m) — > ^^47r 



■ -jrnL _l_ 



(37) 



The double limit A ^ 0, m ^ can also be obtained in the opposite order, and the result is the same: 



6Ao{A,m = 0) = 



i \ - 1^ 



fr2 



-jLko 



- Ci{jLA) sin(jLA) + Si{jLA) cos{jLA) - | cos(jLA) 



gUL,o,A) 



iA^ 



4^ / djjg{jL,0,A) 
Jo 



2^3 



(38) 



The finite-temperature part converges very quickly: 



L3 



E: 



(2^)2 



;2 - A2 



A— 



0. 



(39) 



The second sum we need is 



m) 



PL3 



E 



(2^)4 



(40) 



We use Ea. (|31|l to separate it into a temperature- independent (Bq) and a temperature-dependent part (B^^). The 
first one is with Eq. 132|l 



d*k 



I y-^l f dko 1 

327 7 



(2^)4 (fco -t- iA)2 k^ + P+ ^2 n^L ^ J J 27r (fco + »A)2 7^ fc^ + P + m- 



dk- 



k 



5-e 



■sin(jfcL) , (41) 



5Bo=Bo-B(/3-.oo,L^(X3) = ^^-y^ rffco (fcg + m^)^e-^-^V^^ . (42) 

A'-^h{jL,m.A) 

The sums over j converge, given the asymptotic behaviors 

h{jL,m,0) --Jl^ + .... (43) 

Ea. (|43|l can be obtained from the integral representation above noticing that, for large j, the integral is dominated 
by small values of kg. These relations show that the sum in Ea. (|42l) converges (quickly). 

We also need the value of (5B at A = 0, where h{j, L, m, A) is apparently infrared divergent, but this limit is actually 
finite. It is most easily obtained by continuing the ko integral to 1 -|- e dimensions and taking the e — > limit at the 
end, with Kn again a modified Bessel function: 

oo 

2^2 



1 ^ > 



A 

A7r^L2-.jO{jLm)^J''\:,2_if-'^ 

m 



47^2 _L3 j 



J2 4 [0'^"^ - ULmf)Ko{jLm) + 2(1 + {jLmf)K,{jLm)] , 



We can further take the limit m — > 0: 



27r2L4 ^ ,-4 2L^ ■ 

To take the double limit in the opposite order leads - not surprisingly - to the same result: 

Jt^o V ^ ' 

A3/j(jL,0,A) 

A" 



St 



j¥0 



4 6(jLA)2 _^ 2Ci(jLA) [jLA cos(jLA) + 4sin(jLA)] 



(jLA)3 

+ 2Si{jLA) \jLAsm{jLA) - 4cos(jLA)] 

' 27r2L4 2^ j4 - 2L4- 

After performing the correct analytic continuation, the temperature-dependent part is 

1 ^ u;l + A^ C-j^r 1 
L3^(u;2-A2)2 a;„ e'^'-n-l' 

n 

Finally, we use similar steps for C and D: 

^ ' PL^ ^ (2^no//3)' + (27rn/L)2 + m2 

'fT'li-f-'-' 

/d^k fcfl m2 K2{jmL) 1 a;„ 

(2^0^ fc2 + fc2+„2 + (2;^2^ j2 L^4-e/3-n-l • 



r)^fc2 + fc2+„2 (27r)2i2Z^ j2 
(27rn/L)2 



/3i3 Z-^ (27rno//3)2 + (27rn/L)2 + m2 



7 (27r)4 fc2 + fc2+^2 (2^)2iZ.j^ lOm j+ ^.^^ J + L3 1. ^„ e/^-- - 



For TO ^ 0, one has to be careful with the mode n = 
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APPENDIX C: RELATIVISTIC CALCULATION, REAL TIME FORMALISM AND THE CORRECT 

ANALYTIC CONTINUATION 



In order to verify our procedure to compute the finite temperature corrections of the nucleon mass, we now repeat 
the calculation of the simplest diagram by dispensing of the simplifications due to the use of both the heavy-baryon 
and the imaginary time formalisms, followed by analytic continuation to the real axis. 

The real-time finite temperature formalism (RTF) is another way of (perturbatively) computing finite-temperature 
corrections. As opposed to the more common imaginary-time formalism (ITF), it computes correlators directly in 
real time and continuous frequencies. The Feynman rules are very similar to the ones at zero temperature, except 
that the propagators contain an additional term describing the influence of the thermal medium on the propagation 
of the particles . The pion propagator becomes 



iD{k) = 



iO 



+ 27TnB{\ko\)S{kl - ^ m^), 



(51) 



where nsdfcol) = (e'^l'^ol 

^S{p) = {ip^r + M) 



1) ^ is the bosonic distribution function. The fermion propagator is 

^ 27TnF{\p„\)Sipl-f-M^) 



ko + iO 



pI-T? - A'P + iQ 
2Tinp{\M + ko\)6{kQ), 



(52) 



where npikQ) = (el*^"' + 1)~^ is the Fermi distribution function. The physical origin of these extra terms is the Pauli 
blocking (in the case of fermions) or stimulated emission (in the boson case) caused by the real, on-shell particles 
present in the medium. In the fermionic case, for instance, a state that is fully occupied {np = 1) reverts the sign of 
the "ie" prescription and the fermion can propagate as a hole. Notice that the number density of particles in the 
heavy baryon propagator is np{\M -\- ko\) (as opposed to rii? (|fco|)) and therefore exponentially small at all 
temperatures f3M ^ 1 where the effective theory applies. As an example, let us compute the real part of the second 
diagram in Fig. [3 for notational simplicity in the infinite volume. Up to irrelevant constants. 



iG{E) 



d*k 



fc2((£; + fco)/-fc.7 + MA) 



"■0 



fc2 



tO 



2TTnB{\ko\)S{k^ 



{E 



kaY - fc2 - Ml 



iQ 



+ 2iTTnF{E + ko)6{{E + kof 



Ml) 



Using the relation 



= V 



— iTrS{x), 



(53) 



(54) 



where V stands for the principal value, the real part of G'(w) is 



TZeGiuj) 



dPk 



l + 2nB{uJk) {E + uJk)l°-k-^ + M/s, {E - LJk)l° - k ■ ^ + A4a 



(27r)3 
1 - 2nF{\/^^ 



2ujk 



Ml 



Ml) 



V {E + ujk)^-k^ 
-iJk^ + MIj° - ifc • 7 + A/a 



{E - CJfe)2 - fc2 _ 



E=M+u) 




(55) 



^ In diagrams with more than one loop, the RTF rules are a little more involved. 
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This last result can also be obtained directly using the RTF with the heavy baryon propagator in Ea. H52|) . 

On the other hand, we compute a related quantity, namely the Euclidean time Matsubara function Giuj) defined for 

discrete imaginary values of -E = — Z7r(2n + l)//3, as 



TleGiE) 



d'^k p 

d^k 
(2^ 

d^k 
(2^ 



{E + w„)2 + P + + 

-i{E + kQ)-i'^ -ik-f + Ma 1 



Ml 



2wfc 



i{E + iwfc)7" -ik-f-^ Ma 



-i(£' - «t^fc)7° - zfc • 7^ + Ma 



1 



[E + iCJfc)2 + fc2 + Ml 

P + M|,70 ^ik-j^ + MA 



(E 



2A/fc2 



Mi 



icjfe)^ + fc2 + 
1 



(-£ + ?Yfc2 +M2)2 +cj2 

-iJk^ + MIj" -ik-^^ + Ma 



1 



1 



(E- 



/j2 



M2)2 



Jiy/W+Ml-il3E 



1 



(56) 



We now substitute again e^^^^ by —1. The real part of G{E) is related to the Matsubara function t/(w) by analytic 
continuation: G{E) = —Q{—iE + 0) |^ . The direct calculation of G{E) using the RTF and the indirect one through 
analytic continuation from the ITF agree, and they also reproduce in the limit M ^ oo the calculation using both 
the ITF and the heavy baryon formalism discussed in the main text and in Appendix B. 
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